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1 Preliminary considerations and statement of the results 

Rigorous study of growth processes, and constant attempts to give mathematical 
rigor to Kardar-Parisi-Zhang theory and its predictions, led to a very rich flow of 
results. Remarkable, but yet very limited progress was achieved in the last few 
years for 1 + 1 dimensional models by using a broad spectrum of techniques and 
arguments from the theory of random matrices, first passage percolation and in- 
teracting particle systems. An important role in the successful application and 
interpretation of obtained results was played by the fact that some properties of 

V. Beffara 

UMPA - ENS Lyon, 46 Allee d'ltalie, F-69364 Lyon Cedex 07, France 
E-mail: vbeffara@ens-lyon.fr 

V. Sidoravicius 

CWL Kruislaan 413, 1098SJ Amsterdam, The Netherlands and 

IMPA, Estrada Dona Castorina 1 10, Rio de Janeiro RJ 22460-320, Brasil 

E-mail: v.sidoravicius@cwi.nl, vladas@impa.br 

M. E. Vares 

CBPF, Rua Dr. Xavier Sigaud 150, Rio de Janeiro RJ 22290-180, Brasil 
E-mail: eulalia@cbpf br 



2 



Vincent Beffara et al. 



a particular 1+1 dimensional growth could be mapped and studied in terms of 
the length of a maximal increasing sub-sequence in permutation of A'^ elements, 
with N ^ or as properties of the one-dimensional totally asymmetric simple 
exclusion process, or finally as a directed polymer in two dimensions subject to a 
random potential. 

The starting point and motivation of the present work was the following question: 
how can a localized defect, especially if it is small, affect the macroscopic behav- 
ior of growth system? This is one of the fundamental questions in non-equilibrium 
growth: is the asymptotic shape changed (faceted) in the macroscopic neighbor- 
hood of such a defect at any value of its strength, or, when the defect is too weak, 
then the fluctuations of the bulk evolution become predominant and destroy the 
effects of the obstruction in such a way that its presence becomes macroscopically 
undetectable? 




Fig. 1 Example: effect of the local defect on propagation of slow-combustion fronts: Higher 
concentration of potassium nitrate along vertical line changes shapes of otherwise nearly hor- 
izontal front lines. Higher velocity of propagation of combustion front along this line affects 
macroscopic picture producing a faceted shape around the defected region. (Experiment pro- 
duced by M. Myllys et al. il2|. Reproduced by permission.) 



Such a vanishing presence of the macroscopic effect as a function of the strength 
of obstruction represents what is called dynamic phase transition. The existence 
of such a transition, its scaling properties, the shape of the density profile near the 
obstruction, and also whether information percolates through the obstruction, are 
among the most important issues. As we already mentioned, in many cases these 
questions can be translated into the language of pinning of a directed polymer in 
presence of bulk disorder, or to the question of whether the presence of a static 
obstruction in a driven flow, such as a slow bond in one dimensional totally asym- 
metric simple exclusion process (TASEP), always results in change of the current, 
as predicted by mean-field theory (see |[9l). 
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Similar questions are often raised in the setting of first-passage percolation (FPP) 
(see moil for a review). More specifically, consider an FPP process in which the 
passage times are shortened (in some random or deterministic way) along one 
of the axes. Such type of perturbation usually is referred in physics literature as 
columnar defect. Is it then true that the asymptotic shape of the cluster is always 
modified, or is there a phase transition for the amount of perturbation? Fig. 1 
shows experimentally obtained successive slow-combustion fronts on the piece of 
paper, in the presence of the columnar defect (see 1121 for details). 

All these questions remain rather poorly understood on the experimental and the 
theoretical physics level, and essentially no rigorous mathematical results are avail- 
able Q One particular model, however, known as Poly-Nuclear Growth (or PNG for 
short) makes the exception. In this process, on a flat substrate, layers appear on 
top of each other according to a Poisson process, and then grow laterally with 
constant speed —1 and +1, respectively, and merging whenever they collide. This 
model is interesting in the sense that it lends itself to the explicit computation of 
various quantities, such as growth speed and interface fluctuations (see Il4l ll3l ). 
The proofs depend on exact algebraic computations, and it is not clear how to ap- 
ply them as soon as the dynamics are modified in a non-translation-invariant way 
(with an exception of certain symmetrized cases (see IH)). To treat the PNG case, 
new arguments have been developed in [TJ and 161 . 

The process which we introduce and study in this paper, the Randomized Poly- 
Nuclear Growth (or RPNG for short), has several apparent similarities with PNG, 
but the substantial difference is that after the nucleation, boundaries of layers 
(plateaux) move as independent random walks instead of moving deterministi- 
cally. This dynamics corresponds to the modified Glauber dynamics for the two- 
dimensional Ising model on a solid surface at temperature 0, as considered in [8\, 
with additional nucleations. Informally it can be described as follows: On each 
edge of the discrete line Z, with some positive rate A, pedestals of width and 
height 1 appear. Every pedestal will evolve into a plateau, its endpoints perform- 
ing independent continuous-time, discrete-space symmetric random walks with 
jump rate 1, until one of the following happens: 

- If two plateaux collide, i.e. if the right edge of one occupies the same site as the 
left edge of the next one, then they merge and only the two extremal endpoints 
continue to perform random walks; 

- If one of the two endpoints of a plateau jumps to the location of the other one, 
so that the plateau shrinks to width 0, then it vanishes. 

If a pedestal appears on an edge which is already spanned by one of the 
plateaux, it appears on an additional layer on top of it. Only two plateaux evolving 



' In this context it also interesting to mention a recent work (see 1 2 1), where alternative meth- 
ods to detect localized inhomogeneity were used. In |2| the following question was studied: 
Consider a graph with a set of vertices and oriented edges connecting pairs of vertices. Each 
vertex is associated with a random variable and these are assumed to be independent. In this set- 
ting, the authors consider the following hypothesis testing problem: under the null, the random 
variables have common distribution A'{0, 1), while under the alternative, there is an unknown 
path along which random variables have distribution N{ii, 1), /x > 0, and distribution A'(0, 1) 
away from it. The authors study for which values of the mean shift fi one can reliably detect the 
presence of such defected path and for which values this is impossible. 
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within the same layer can merge, and besides the growth of the plateaux at differ- 
ent layers has to obey the following constraint: Every plateau which is not on the 
ground level must always lie entirely on top of a plateau at the layer below it. Each 
jump which would break this constraint is simply suppressed. For any fixed value 
of A (see Proposition [T] below) the process has linear growth, i.e. if by ht{e) we 
denote the height of the aggregate at time t above the edge e, then the following 
limit exists: 

v{k):= \\mJy{e)/t>Q. (1) 

Notice that in the case of a strictly growing process (such as first- or last-passage 
percolation), the law of large numbers is usually a direct consequence of a state- 
ment of sub-additivity, while if the dynamics is allowed to both grow and shrink, 
like in the present case, even the existence of the asymptotic speed is non-trivial. 

Next we will introduce a "columnar defect" to the system by changing the 
nucleation rate on the edge eo :=< 0, 1 > from A to A + Ao, with Ao > 0, and the 
main question which we address is how does this affect the asymptotic growth 

rate of the aggregate above the origin. More precisely, let /i, (eo) be the height 
of the aggregate above eo at time t under the modified dynamics (we make the 
dependency on A implicit in this notation, since the parameter of interest is Ao). 
Assuming that the following limit exists 

v(A,Ao):=^lim^/7^(eo)A, (2) 

we are interested in whether the inequality 

v(A,Ao)>v(A) = v(A,0) (3) 

holds for all Ao > and set: Ac(A) := inf{Ao : v(A,Ao) > v(A)} be the critical 
value of Afl. 

We are now ready to state our main result. 

Theorem 1 

v(A,Ao) := ^lirn^^^!^ = A +max{Ao - 1,0} 

almost surely. In particular, Ac(A) = 1, for any value of Pi. 

In proving this theorem, we actually obtain a more precise description of the 
behavior of the system in the supercritical phase Ao > Ac (A): the behavior of 

h^{e) for e 7^ eo is perturbed only sub-linearly. In other words, the defect cre- 
ates a kind of "antenna" at the origin. This contrasts with the expected behavior of 
the usual PNG with a supercritical localized defect, for which the region of influ- 
ence of the additional nucleation process is believed to extend linearly in space. 
However if in the RPNG model we would remove assumption of the symmetry of 
walks and would take the boundaries of plateaux moving as random walks with 
a drift (left boundary with a negative drift and right boundary with positive one), 
we expect the behavior of the system to resemble the behavior of PNG. 

Moreover, the phase transition is non-trivial even without nucleation in the 
bulk (i.e. A = 0). Even in this simplified situation, for any value Ao > the origin 
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is eventually covered by plateaux for all times, however < := Aj.(0) < +°° so 
that if Ao < then v(0, Ao) =0, i.e. there is sub-linear growth, while for Ao > A^ 
we have v(0, Ao) > 0. As we shall see. A/, = 1, in agreement with Theorem 1. 

Comment It might be instructive to compare with the following simplified model 
in the half-line ]R+. Consider a situation corresponding to A = and an additional 
vertical wall at vertex to which the left endpoint of each pedestal is fixed: the 
initial configuration is empty; at each occurrence of a Poisson process in M+ of 
intensity Ao, a new pedestal is created at the edge eo — and this is the only source 
in the system. Right boundaries of pedestals move according to the same mecha- 
nism as described in the basic homogeneous model, except that now each pedestal 
is annihilated upon its right boundary hitting the site 0. 

Of particular interest is the particular case Ao = 1 /2, when the system is in ex- 
act bijection with 0-temperature Glauber dynamics for the two-dimensional Ising 
model, with an initial configuration consisting in the first quadrant uniformly at 
spin — 1 and the rest of the lattice uniformly at spin + 1 . The height function at eo, 
here represented by Nt, the number of pedestals in the system at time t, describes 
the shape of the +/— interface, and in the large time limit it is known to behave 
according to Lifshits' law of motion by mean-curvature (see ID, |Il4 |). 

In this particular case, there is also a one-to-one correspondence with the sim- 
ple symmetric exclusion process: A', corresponds to the position of the rightmost 
particle in a symmetric simple exclusion process in Z starting with the configura- 
tion where the occupied sites are exactly those x < 0. A classical result of Arratia 
(see Theorem 2 in |3|, cf. also [11], comments at p. 414) tells in particular that 
A^;/f — > almost surely. A simple comparison with this case tells us that 1/2 is the 
critical value of Aq in this simplified model. 

Structure of the paper. In Section|2]we give a formal construction of the process, 
in terms of a Harris-type representation, and discuss two other equivalent repre- 
sentations. In Section [3] we derive basic properties of the homogeneous process 
without the defect at eo- More specifically, we prove a law of large numbers for 
the height of the interface above edge eo; a key step towards this result is the proof 
of eventual covering of any finite-size interval of the line. Finally, in Section|4]we 
present the proof of Theorem[Tl 

2 Construction of the process: three ways to look at the system 

2.1 Description as an interacting particle system 

One can describe this process in terms of an interacting particle system. Each left 
endpoint of a plateau is represented by a left-particle or ^-particle for short, at the 
corresponding site of Z, and each right endpoint of a plateau as a right-particle, or 
r-particle (which can be interpreted as an anti-particle). If several plateaux have 
their endpoints at the same site, this site will contain several particles, but notice 
that a site will never contain particles of both types. 

The dynamics then looks like a zero-range process with annihilation and nu- 
cleation, i.e.: 

- For each non-empty site, with rate 1, one of the particles jumps with equal 
probability to a neighboring site; and if the target site contains a particle of the 
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Other type, then the jumping particle and one of the target particles annihilate 
instantaneously; 

- On every edge (x,x+ 1) of Z, with rate A, there is nucleation, meaning that 
an ^-particle appears at x and an r-particle appears at x + 1; in each case, the 
particle is instantaneously annihilated (together with a particle of the opposite 
type), if a particle of the opposite type is already present at the corresponding 
site. 

Remark. Notice that the particle configuration at time t gives only partial infor- 
mation about the initial system; in particular, it gives no way of following the 
evolution of a given plateau boundary as a function of time — indeed, with our 
description of the particle dynamics, the notion of a tagged particle is ill-defined. 
Nevertheless, knowing the particle configuration at time t allows for the recon- 
struction of the plateaux configuration at the same instant up to a global height 
change; and of course, the full space-time history contains exactly the same infor- 
mation for both descriptions of the system. 

We first construct the above described interacting particle system using Harris 

graphical representation. Let E : = {< 1 >,x e Z} be the set of edges of the 
integer line Z. Define the following three independent Poisson systems: 

- £^r = {^r{e),e e E} is a collection of independent Poisson processes of left- 
arrows, each with intensity dt/2; the individual arrows (or marks) of the pro- 
cess at edge e will be denoted by ^"(e), n ^ 1, with similar notations for the 
other two systems; 

- ^1 = {^f(e),e e E} is a collection of independent Poisson processes of left- 
arrows, each with intensity dt/2; 

- ^r = {^r (e) , e € E} is a collection of independent Poisson processes of right- 
arrows, each with intensity dt/2; 

- ^1 = {Mt{e),e e E} is a collection of independent Poisson processes of right- 
arrows, each with intensity dt /2; 

- ^ = {^,eeE}isa collection of independent Poisson processes of undi- 
rected edges, each with intensity Xdt. 

££r (respectively will correspond to jiunps of right particles (respectively left 
particles) to the left, Sir and to jumps of right and, respectively, left particles 
to the right, and jV to nucleation times. 

Given 7" > 0, we first construct the process up to time T , as follows. There 
is positive density of edges carrying no mark of any of the Poisson before time 
T; between any two such edges, there are almost surely only finitely many marks 
appearing by time T, so that they can be ordered chronologically, and we only 
need to describe the action of each of the marks. 

- If there is a directed arrow of type r or, respectively, type I at edge e and time t, 
and if at this time the site at the source of the arrow is not empty and contains 
at least one right or, respectively, at least one left particle, then one particle 
at this site jumps in the direction of the arrow; moreover, if after this jump it 
meets a particle of the opposite type, they annihilate instantaneously. 

- If there is a mark of J/ at edge e=<x,x+l > and time t, then at this time 
an ^-particle appears at site x and an r-particle appears at site jc -I- 1; in each 
case if the corresponding site contain particles of the opposite type, a pair of 
opposite type particles is annihilated instantaneously. 
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It is then easy to see that the obtained processes for different values of T are 
consistent, and letting T go to infinity gives the process for all time. We will refer 
to this version of the construction as the process in quenched environment, i.e. 
seeing the Poisson processes as a random scenery in which a deterministic process 

evolves. 

Remark. Just for the sake of constructing of the process, there is no need to distin- 
guish between and J^^, or and The distinction plays a role only when 
describing the evolution of the symmetric difference between two configurations, 
which we shall need later. 

Particle system corresponding to the process witli a colmnnar defect: A colum- 
nar defect in the plateaux dynamics, discussed in the Introduction, is described in 
the interacting particle system formalism by the addition an auxiliary nucleation 
process J^o with rate Ao > localized on the edge eo- In other words, the pedestals 
appear with rate X everywhere except at eo where the nucleation rate is A + Ao. 
The construction of the new process can be done in exactly the same way as for 
the unperturbed system, using a graphical construction and letting the marks of 
^ act the same way as those of 

Notation: nx(t) will denote the number of particles at site jc and time t in the above 
construction; .jVe{t) denotes the number of nucleation marks occurring at the edge 
e up to time t. 

2.2 Construction in terms of interacting random walks 

In this subsection we will give an alternative interacting-random-walks-type de- 
scription of our system. To do so, one needs a rule to determine how do these 
walks interact, and which of the particles at a given site will jump when an arrow 
(in the previous construction) appears, and which will be annihilated when a parti- 
cle of the opposite type arises. Obviously the process itself will not be affected by 
the choice, so we have some freedom here. Some of the more natural choices are 
to declare that the oldest particle {i.e. the one that appeared earliest in the system), 
or the youngest particle, or the one that has been at the given site for the longest 
time, or the shortest time, is "active". The rule can even be different for jumps and 
annihilations. We will refer to such a rule as a priority rule. 
In particular, the previous graphical construction can be replaced with the follow- 
ing interpretation of the dynamics: 

1) Evolution of a single particle. Whenever a particle i is bom into the system, 
it comes together with a "project", i.e. a family (7^") of independent exponential 
times of rate 1 and a family (e") of independent signs, and a clock. Those will 
represent waiting times between consecutive jumps, and, respectively, directions 
of jumps for the path of this particle, if it were moving alone in an otherwise empty 
space. 

2) Interaction. Choose any priority rule. When a particle is alone, it is always "ac- 
tive", i.e. its clock runs deterministically with speed 1, and it executes its project, 
i.e. jumps according to times and directions prescribed in its project. When there 
are several particles at a given site, only one of them is active, according to the 
chosen rule, and the others are "sleeping". The clock of the active particle runs 
deterministically with speed 1, the clocks of the other particles do not run at all. 
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Whenever the clock of any particle reaches the next waiting time for the project 
of this particle, the particle jumps and the clock is reset. 

Remark. Another natural choice for a jumping rule could be that at a site with n 
particles, the clock of each of them runs with speed 1/n; again the distribution of 
the numbers of particles at different sites would be the same as in the graphical 
construction. 

The formal construction of this dynamics could be done similarly as in the Harris 
construction: Given a space-time configuration of particles (space time Poisson 
process, corresponding to the .A' system in the graphical construction) and their 
families of "projects", for any fixed T > Q, there are infinitely many edges to the 
left and to the right of the origin such that by time T they will not be crossed by any 
type of particle, neither a nucleation will occur at them by time T . It follows quite 
simply that the evolution of the system within the space time block in between two 
such edges and [0, T] depends only on finitely many particles, being independent 
from the evolution outside. 



2.3 Back to the PNG-like description 

We saw how to go from a PNG-like description to a family of interacting parti- 
cles, but our final goal is to describe the long term behavior of the height {i.e. the 
number of layers) at a given edge, so we need a way to go back to the layer struc- 
ture. As we saw already the f-particles correspond to left boundaries of plateaux, 
and the r-particles correspond to right boundaries. It only remains to compute the 
height at edge eo =< 0, 1 >. 

Let (^/) be the flow of ^-particles through the edge e^. (/>o = and 0^ increases 
(resp. decreases) by 1 whenever an ^-particle crosses eo to the right (resp. to the 
left). Define ((/)/') similarly as the flow of r-particles. Let /i,(eo) be the height of 
the profile at edge eo and time t+ {i.e., the cadlag version in time of the height 
process). Then, recalling previous notation, 

/jr(eo)=^o(0 -</>/ + </>;, (4) 

From ht (eo) and the previous discussion, it is then easy to obtain the whole profile 
of plateaux. 

Remark. The graphical method gives us a construction of h, := {/2,(e)}ggE as a 
Markov process. The {/j,}-evolution is attractive (or monotone in the sense of 
definition 2.3 of Chap. II in (TT\), though the particle system is not. The height 
process is also monotone in the nucleation rates. 



3 Behavior of the homogeneous process (A > 0, = 0) 

3.1 Level lines, Hammersley-process-like description 

Let L, (resp. Rt) be the position of the left (resp. right) end of the first-layer plateau 
containing eo at time f, if there is one i.e. if ht{eo) ^1 — for completeness, one 
can let Lt = +°° and Rt = ~°° if h,{eo) = 0. 
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Let T be the first time at which R, = Lf + 2, if any (otherwise T = +°°). In the 
absence of other layers and of further nucleation, the behavior after T of L, and 
R, until the first time when 7?, = L, + 1 would exactly be that of two independent, 
centered continuous-time random walks. The presence of the second and above 
layers, by the previous description, can only affect the pair by the suppression 
of left jumps of Rf and right jumps of L, — notice that if, in the description at 
Subsection 2.1, there are several particles at R, (L,) at time t and one has to jump 
to the right (left, resp.), then so does R, {L, resp.). 

Moreover, with rate X there is nucleation on the edge < 7?, , 7?, + 1 >, and when 
this occurs, 7?, jumps to the right by at least one step. The last thing that can affect 
7?; is if the plateau [Lf,7?,] merges to the next plateau to its right, which again can 
only increase the value of Rf. The same remarks apply mutatis mutandis to Lf. 

Taking all those effects into account, one obtains that the process (7?,) stochas- 
tically dominates the biased random walk {Rt) that has rates 1 /2 to jump to the 
left and 1/2 + A to jump to the right. Similarly, {L,) is stochastically dominated 
by [Lf) which is biased to the left. Since the probability 



is positive, using standard arguments we see that for to large enough we have 



for suitable positive constants c\,C2- Applying the Borel-Cantelli Lemma, one in 
particular sees that a.s. any compact interval of the real line will eventually remain 
covered by the plateau [Lt , 7?(] . Moreover, extending the same line of reasoning one 
sees that 



for suitable constant C > 0. 

We shall give an additional description of the time evolution of the system, in 
terms of a family of space-time curves similar to a Hammersley process (cf. HI). 
The main difference is that the paths of the particles between meeting points are 
random walks instead of straight lines. 

More precisely, one can see the path of any given particle, from its birth to 
its eventual annihilation, as a sequence of vertical and horizontal line segments. 
Let be the union of the space-time paths of all the particles corresponding to 
endpoints of plateaux of the first layer, together with the segments joining the birth 
places in space-time of the corresponding pairs. 

As we have seen above, for each bounded segment [— A',^r] there is an a.s. fi- 
nite random time tK, such that, the segment [—A', K] will be covered by a plateau 
of the first level for all t >tx- This means that is the disjoint union of a count- 
able family of finite loops (corresponding to plateaux with a finite life span) and 
one unbounded curve, which we will denote by Ci and call the first level line of 
the process. 

Similarly, let 'i^j^ be the union of all the space-time paths of the endpoints of 
plateaux of the ^th layer together with the segments joining the birth places in 



p=P{yt>0, R,-L, > 1) 



(5) 




(6) 



liminf — log7?, > C a.s. 
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Fig. 2 Illustration of an RPNG process. (Level lines are bold.) 

space-time of the coiresponding pairs, "^j. is then the disjoint union of a count- 
able family of finite loops, and one unbounded curve Q. Notice that in general a 
component of "^j. and a component of % for I need not be disjoint. 

Using this description, it is possible to compute the height of the RPNG pro- 
cess at any edge e and (almost) any time t: Indeed, it is equal to the number of level 
lines that separate from the line f = {i.e. that pass below it an odd number 
of times), plus the number of loops containing (e,f). This means that the law a 
of large numbers for h, can be translated into a statement about the asymptotic 
density of the level lines. 

3.2 The asymptotic speed 

The first step in proving a strong law of large numbers for the homogeneous pro- 
cess is to obtain an a priori bound for the number of particles present at a given 
site, at a given time. We obtain such a bound in the form of a stochastic domination 
result. 

Lemma 1 There exist three countable families {X^)xsZ, {Xr,.)x€Z ond {Y!')xs'L of 
integer-valued stochastic processes, and a coupling of them with the homogeneous 
particle system introduced above satisfying the following conditions: 

1. For each x G Z, {X-^^)i^q and {X^,)t^o are simple symmetric random walks on 
Z+, with reflection at the origin, jump rate 1 at the origin and 2 elswhere, and 
(y,-*)(j,o is a simple symmetric random walk on Z with jump rate X; 

2. Almost surely, for all x Gl^ and t ^0, 

nx{t)^Xl,+X^,+ sup {Y:~Y^,). (7) 

s,i'6[0/] 

Proof Let jc G Z, and in the graphical construction of the process (cf. Section |2]l, 
let {T( x.k)k>o (resp. {Tf^f,)) be the kth time at which there is an £-arrow pointing 
toward (resp. away from) site x. Define X^ as follows: XJ^ = 0, X^ jumps up by 1 
at each Tix,k, and jumps down by 1 at each t G {Tf^/.} such thatX^j_ > 0. Clearly, 
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defines a continuous-time random walk as in the statement. Similarly we define 
using r-arrows. Let Y/' = ^i^x.x+\){l) ^ -^.c- 1 .a) (0 • 

rtx can change for two reasons (nucleation or jumps): when a particle at x 
jumps away, % decreases by 1, and when a particle at x± 1 jumps onto x, it ei- 
ther increases by 1, or (if Hx > 0) it decreases by 1 depending on the types of 
particles involved; if there is nucleation at the edge {x — l,x) or {x,x + I), then % 
again increases or decreases by one according to the types of particles involved, 
as before. 

Fix t >0, and let 

T := sup{5' < t : nx{s) = 0}. 

If T = f and nx{t) = we have nothing to prove. From now on we are assuming that 
this is not the case. In particular, during the time interval [t,?], the site x is never 
empty and the type of particles present in it does not change; besides, nx{T—) = 0. 

Apply the particle-counting of the previous paragraph between times t and t. 
The contribution from nucleations is exactly equal to ±(F/' — Y!^_) (where the sign 
depends on the type of particles present at site x), and the contribution from jumps 
is controlled by the jumps of Xf^ or X^^. This leads to the following estimate: 

nx{t)^\Y?-Y',-\+Xl+X',,. 
The bound in the statement of the Lemma directly follows. 
Proposition 1 Recall that cq =< 0, 1 >. For every A > 0, almost surely, 

limMf^l^A. 

t^=<= t 

Proof Let Q > be a positive integer, and let 

HtiQ):^ ht{<x,x+l>) 

x=-Q 

be the integrated height of the process over the interval [—Q, Q\. The variations of 
Ht{Q) in time come from several sources (according to the definition of the height 

process): 

- Whenever there is a nucleation on one of the edges in [— 2, 2] , Hf {Q) increases 
by 1; 

- Whenever an ^-particle in[— Q+1,Q— 1] jumps to the right (resp. to the left), 
Ht{Q) decreases (resp. increases) by 1 

- Whenever an r-particle in [—Q + 1 , g — 1] jimips to the left (resp. to the right), 
Ht{Q) decreases (resp. increases) by 1; 

- Whenever a particle sitting at 2 (resp. —0 jumps to site Q— 1 (resp. —Q+V), 
H({Q) changes by -|-1 or — 1 according to the type of particle (i or r) involved. 

Any other jump or nucleation might change the value of n±Q(f), but it will never 
affect the value of Notice in the second and third cases, the decrease of Ht{Q) 
occurs independently of the jtimp leading to an annihilation or not. 

Notice that a particle sitting in [—2+1,2—1] jumps to the left and to the 
right with the same rate (be it an £- or an r-particle), and that the total jump rate 
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of particles in this interval is, by construction, bounded by 2Q — 1 (the total jump 
rate of particles at any given site being 1 if the site is non-empty, and of course 
if the site is empty). Hence, the contribution of such jumps to H,{Q) is at most of 
the order of the position at time f of a symmetric, continuous-time random walk 
jumping with rate 2Q — 1 — i.e. , they are of order \/t. 

What remains is the contribution of nucleations inside [— Q, which have 
total rate 2QX and satisfy a strong law of large numbers, and the fourth term 
which, although difficult to estimate precisely due to its dependency on the exact 
particle distribution, only contributes to Hi{Q) with rate at most 1 {i.e. , the sum 
of the jump rates from —Q\.o—Q+\ and from QioQ 1). 

Summing all the above contributions to Ht{Q), and then letting t go to +oo, we 
obtain with probability 1 ; 

2eA-l ^liminf^^ sClimsup^^ <2eA + l. (8) 

t ,^oo t 

Now, we know that at time t, the difference in height between two neighboring 
edges, say < x — l,x > and < x,x+l >, only depends on the number (and type) 
of particles present at site x at time f . By (|7]l we see that almost surely, for every 

xez, 

t-.oo t 

The height ht{e) at any edge e G [— (2, Q] can be computed from the height 
at edge eo by only looking at the values of the for x E [—Q+ 1]- In 

particular, for such an e, 

G-l 

\h,{e)~h,{eo)\^ £ K(0|. 

x=-Q+l 

Summing over e to obtain Ht{Q), and taking a crude upper bound, leads to: 

Q-l 

mQ)-2Qhr{eo)\^2Q £ |n,(f)|. 

x=-Q+l 

Combining this with equation ^ and the previous estimate on the n^{t), we ob- 
tain, still with probability 1: 

1^1- • r2e/f (eo) , ,. 2e/j,(eo) . 
2Qa — 1 < limmf < hmsup < 2Qa + 1 

^ — t t — *oo i 

(notice that we are still working with fixed Q). Dividing by 2g, we get 

A < limmf — ^ — - < limsup — ^ — - < A H . 

2Q^ t " t^o.^ t " 2Q 

Since this holds for every positive Q, the proposition is proven. 
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4 Proof of Theorem 1 

We start with a homogeneous RPNG process as described in Section |2l let ht be 
its height function. Notice that the whole dynamics can be described in terms of 
this function. Notation: e-f = {x,x+ 1), for each integer x. 

Now add an extra nucleation with rate at the edge eo to the Harris graph 
of Section 2, and let h'^ be the height function of the modified process, with 
being the number of particles at x in the corresponding new particle system. By 
monotonicity of the plateaux dynamics (cf. Remark right after equation we 
have that h', ^ ht; and let h, = — ht ^ be the discrepancy process. We are 
interested in differences in the law of large numbers between the two versions of 
the process, i.e. in the asymptotic properties of ht jt. 

The discrepancy process can be thought as a particle system, composed of 
what we call virtual particles: According to the difference 

ht{.ex-\)-ht{ex) 
being positive or negative we say that at time t there are 

virtual right or, respectively, virtual left particles at site x. In this interpretation the 
creation of virtual particles corresponds to a nucleation in jV^^ \ jV^^ . 

Using the Harris description of subsection 2.1 (or the random walk represen- 
tation of subsection 2.2) one can check explicitly that the spatial displacements of 
virtual particles remain as of independent random walks even in presence of the 
interaction with other virtual particles while jump times can be possibly delayed 
— or in other words, that conditionally on the whole particle system, the current 
rates of jump of a given virtual particle to the left or to the right always remain 
equal. Doing so in full detail involves looking at many different cases, but the 
discussion can be shortened considerably by the following considerations: 

- If a site contains only virtual particles, then it is obvious from the construction 
that the first jump of one of them away from the site occurs with equal rate to 
the left and to the right; so we only need to check what happens when normal 
and virtual particles share a site. 

- If a site contains both normal and virtual ^-particles, then the perturbed system 
also contains ^-particles. In particular, the original and perturbed system will 
use the same arrows in the Harris construction, meaning that their difference 
will not be affected by the dynamics and the virtual particles will not move as 
long as there are normal particles present, at which point we are back to the 
previous case. Of course, this also applies to the case of r-particles. 

- If a site X contains a > normal ^-particles and > virtual r-particles, then 
one needs to check 3 sub-cases (a < b, a = b and a > b), leading to different 
behavior in the perturbed system. We treat one of them fully, the others being 
similar: Assume that b > a (see Fig. |3]l, so that the perturbed system contains 
{b — a) r-particles at site x. The occurrence of an arrow from x affecting left- 
particles will affect the unperturbed system and not the perturbed one, meaning 
that one of the virtual particles will follow the jump as well; the occurrence of 
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an arrow affecting right-particles will only modify the perturbed system, and 
again in this case a virtual particle will follow; in both cases, the symmetry of 
the arrow processes translates into symmetry for the jumps of virtual particles. 



I 



^ ^ r 

Fig. 3 Comparison of tiie perturbed and unperturbed dynamics in the case h > a.In plain lines 
are the height functions for the unperturbed system, and in dotted lines are the height functions 
for the perturbed one. On top is the initial configuration, and below are configurations after 
occurrence of an arrow in -^r, and J^, respectively. The circle on the axis indicates the 
location of site x. 



Proof of the Theorem. Due to the monotonicity, the height of the process with ad- 
ditional nucleation at the origin stochastically dominates that of the homogeneous 
process. In particular, we already know that almost surely 

liminf^^A. 

Step 1. Define 

A = /J^ (eo) - max{/j; (e_ 1 ) , /jj (e 1 ) } . 

If Ao > 3 + A, then 

liminf — > 0. (9) 

Indeed, if > 0, the following can happen: 

- if h'i{e-i) =h'f{ei), then At may increase by one only because of a nucleation 
at the edge eo, which happens at rate A + Aq, and it decreases by one either 
when some particle jumps away from vertex or 1, which happens with rate 
2, or when nucleation occurs at one of edges e_i or ei, which happens with 
rate 2A, or if a particle jumps from vertex —1 (resp. from vertex 2) to (resp. 
to 1), which happens at rate 1 (if such a particle is present); 

- if h'l (e_ i) <h',{e\), then A, may increase by one either because of a nucleation 
at the edge cq, which happens at rate A -I- Ao, or because of a jump of a particle 
from vertex 2 to vertex 1, provided a particle exists at vertex 2, and which 
happens at rate 1/2. On the other hand A, decreases by one either if a particle 
jumps away from vertex 1, which happens at rate 1, or if a particle jumps from 
vertex to 1, which happens at rate 1/2, or, finally, if nucleation occurs at the 
edge ei, which happens at rate A. 
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Examining similarly the cases when A, =0 and A, < 0, we obtain that At increases 
by one with rate at least A + Ao and decreases by one with rate at most 3 + 2A, and 
(|9]l easily follows. 

Remark. In what follows we want to consider the unperturbed system as a random 
environment on which the virtual particles will evolve. 

Step 2. From Proposition [T] and ^ it follows that if Ao > 3 + A, then 

liminf-(A,(eo) — niax{^,(e_i)^,(ei)}) > 

t—><x> t 

and, in particular there exists < a.s. such that the number «o('*) of left 
virtual particles at the vertex and the number n\{s) of right virtual particles at 
vertex 1 is positive, for all s > T^. This also means that no virtual left (resp. right) 
particle created after enters [2, +oo) (resp. (-oo, -l]). 

Step 3. If Afl > 3 + A, then the total number Y,k=2'^ki^) °f virtual (in particular 
virtual right) particles in [2, +oo) grows sub-linearly in f. Indeed, if Ao > 3 + A, 
then from Step 2 it follows that there exists < +oo a.s. such that 4j > for all 
s >T^. After time the system is essentially divided into two half space systems. 
A virtual right-particle born after time (necessarily at site 1) will then follow 
the path of a simple random walk, possibly slowed down by the environment, until 
it first reaches site 0, when it will be annihilated; in particular, hitting probabilities 
for virtual particles are given by a standard gambler's ruin estimate. 

Let A" be a positive integer. Consider all virtual particles in [2, +oo) at large 
times t >Tif. They can be subdivided into three classes: 

- Those which were already present at time T^,, of which there are finitely many; 

- Those which were created after time and entered the half-line [K,+o°) 
by time t. The number of such particles is almost surely asymptotically in t 
bounded above by (1 +e)t/K, since the environment can possibly slow vir- 
tual particles down, but does not affect their trajectories. 

- Those which were created after time T*, but never exited the box [1,A']. From 
Lemma 1 we get that the number of particles at any given site (except for 
and 1) behaves sub-linearly — notice that, since the proof of Lemma 1 is 
purely local, it still applies here to any site that is not directly affected by the 
additional nucleation at eo, which is the case here. Thus, the number of such 
particles at any x in the box is a.s. less than t /K^ for all t large enough. In 
particular, the total number of such particles is smaller than t/K for all t large 
enough. 

Thus, taking K arbitrarily large we obtain the statement. 

Step 4. Observe, that for any t < there exists < R{t) < +°°, random, such 
that h,{< x,x+ 1 >) = if \x\ > R. Thus, if Ao > 3 + A, Step 3 implies that ht{e) 
grows sub-linearly in time for all e 7^ eo- Putting this together with Proposition [T] 
we get that 

lim — J — =A, 67^60, 

for such values of Ao. By the monotonicity of the height dynamics, this immedi- 
ately extends to all Ao > 0. 
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Step 5. Applying the same line of reasoning as in Step 1 we see that max{At , 0} 
grows at most sub-linearly if Ao < 1. Putting together with Step 4, it follows that 

lim^^A, ifAo<l. 

step 6. The height process h[ (eo) grows by one with rate at least A + Aq, and can 
only decrease if a particle (of the appropriate type) jumps across the edge eo — 
which happens with rate at most 1 . This implies that, if Aq > 1 , it is bounded below 
by a random walk with drift A + Ao — 1 . On the other hand, the conclusion of Step 
4 implies that for all e ^ cq, h[{e) satisfies a law of large numbers with speed 
A < A + Afl — 1, thus for all t large enough n',{0) and nj(l) are strictly positive, 
and so 4, as well. Putting these two facts together and using the same reasoning as 
above, we see that the increase and decrease rates for h[ (eo) become exactly A + Ao 
and 1, respectively, for large t, which implies the stated law for large nvmibers for 
A;(eo),ifAo> 1. 

This ends the proof of the main Theorem. 
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